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ABSTRACT 

We discuss a series of observations of the black hole candidate GX 339-4 in low luminosity, spectrally hard 
states. We present spectral analysis of three separate archival Advanced Satellite for Cosmology and Astrophysics 
(ASCA) data sets and eight separate Rossi X-ray Timing Explorer (RXTE) data sets. Three of the RXTE obser- 
vations were strictly simultaneous with 843 MHz and 8.3-9. 1 GHz radio observations. All of these observations 
have (3-9 keV) flux ;$ 10 -9 ergs s" 1 cm -2 . The ASCA data show evidence for an « 6.4 keV Fe line with equiv- 
alent width « 40 eV, as well as evidence for a soft excess that is well-modeled by a power law plus a multicolor 
blackbody spectrum with peak temperature « 1 50-200 eV. The RXTE data sets also show evidence of an Fe line 
with equivalent widths to 20-100eV. Reflection models show a hardening of the RXTE spectra with decreasing 
X-ray flux; however, these models do not exhibit evidence of a correlation between the photon index of the in- 
cident power law flux and the solid angle subtended by the reflector. ‘Sphere+disk* Comptonization models and 
Advection Dominated Accretion Flow (ADAF) models also provide reasonable descriptions of the RXTE data. 
The former models yield coronal temperatures in the range 20-50 keV and optical depths of r % 3. The model fits 
to the X-ray data, however, do not simultaneously explain the observed radio properties. The most likely source 
of the radio flux is synchrotron emission from an extended outflow of extent greater than C?(10 7 GM/c 2 ). 

Subject headings: accretion — black hole physics — Stars: binaries — X-rays: Stars 


1. INTRODUCTION 

The galactic black hole candidate (BHC) GX 339^4 is 
unique among persistent sources in that it shows a wide vari- 
ety of spectral states and transitions among these states. In pre- 
sumed order of increasing bolometric luminosity, GX 339-4 
exhibits states with hard, power-law spectra (‘off state*, Ilo- 
vaisky etal. 1986; ‘low state’, Grebenevetal. 1991); a soft state 
with no evidence of a power law tail (‘high state*; Grebenev 
et al. 1991); and a very bright, soft state with extended power- 
law tail (‘very high state*; Miyamoto et al. 1991). There also 
are apparently times when the flux is high, but the spectrum is 
not as soft as the ‘high* or ‘very high state*. M6ndez & van 
der Klis (1997) refer to this as the ‘intermediate state*. We also 
note that there is some evidence of overlap between the states. 
The broad-band (GRANAT/SIGMA) hard state data presented 
by Grebenev et al. (1991) apparently represents a more lumi- 
nous state than does the broad-band soft state data taken with 
the same instrument. [Miyamoto et al. (1995) has suggested the 
possibility of hysteresis in galactic BHC state transitions.] Sim- 
ilarly diverse sets of states have been observed in X-ray tran- 
sients such as Nova Muscae (Kitamoto et al. 1992; Miyamoto 
et al. 1994); however, GX 339^4 is closer to being a persistent 
source. 

Although there have been a number of observations of 
GX 339-4 in the near-infrared and optical (Doxsey et al. 1979; 
Motch et al. 1983; Motch et al. 1985; Steiman-Cameron et al. 
1990; Imamura et al. 1990; Cowley et al. 1991), including de- 
tection of a 14.8 hr periodicity in the optical (Callanan et al. 


1991), there is no convincing mass function for the system. In 
the optical, the system is faint, variable (My to 1 6-20), and red- 
dened (Ay = 3.5). The physical source of the optical emission 
is unknown. It has been hypothesized that it is entirely dom- 
inated by the accretion disk, as the optical flux is apparently 
anticorrelated with the soft X-ray emission (Steiman-Cameron 
et al. 1990; Imamura et al. 1990). These properties of the emis- 
sion have made it difficult to obtain a good distance measure- 
ment, with estimates ranging from 1 .3 kpc (Predehl et al. 1991 ) 
to 8 kpc (Grindlay 1979), with many authors choosing 4 kpc 
(Doxsey et al. 1979; Cowley, Crampton & Hutchings 1987; 
Zdziarski et al. 1998). A careful study of these distance esti- 
mates is presented by Zdziarski et al. (1998). 

GX 339-4 also has been detected in the radio (Sood & 
Campell-Wilson 1994), and possibly even has exhibited ex- 
tended emission (Fender et al. 1997). Within the hard state, the 
radio spectrum is flat/inverted with a spectral index of a = 0. 1- 
0.2 (Fender et al. 1997; Corbel et al. 1997), where the radio 
flux density S v oc i/°. Furthermore, in this state the radio flux is 
correlated with the X-ray and gamma-ray fluxes (Hannikainnen 
et al. 1998), but the radio flux disappears as GX 339-4 transits 
to a higher X-ray flux/softer state (Fender 1998, priv. comm.), 
comparably to Cyg X-l (Pooley, Fender & Brocksopp 1998). 

During the Rossi X-ray Timing Explorer (RXTE) Cycle 2 
observing phase (1996 December - 1998 February), we per- 
formed a series of eight RXTE observations of GX 339-4. The 
first three observations were spaced a week apart from one 


1 In$titut fur Astronomie und Astrophysik, Abt. Astronomic, Waldhfiuser Str. 64, D-72076 Tubingen, Germany; wilms@astro.uni-tuebingcn.de 

2 JILA, University of Colorado, Campus Box 440, Boulder, CO 80309-0440, USA; mnowak@iocinante.colonido.edu 

^current address, CASA, University of Colorado, Campus 389, Boulder, CO 80309-0389, USA; dove@casa.colorado.edu 

4 Astronomical Institute ‘Anton Fannekoek’, University of Amsterdam, and Center for High Energy Astrophysics, Kruislaan 403, 1098 SJ, Amsterdam, The 
Netherlands; rpf@astro.uva.nl 

^Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK; tiziana@ast.cam.ac.uk 

6 AXAF Fellow; current address Center for Astrophysics, Cambridge, MA 01238, USA; tdimatteo@cfa.harvard.edu 

1 



2 


Low Luminosity States of the Black Hole Car didate GX 339-4. 1. 


another from 1997 February 4 to 1997 February 18. These 
three observations were scheduled to be simultaneous with 8.3- 
9.1GHz radio observations that were conducted at the Aus- 
tralian Telescope Compact Array (ATCA). The results of the 
radio observations have been reported by Corbel et al. (1997). 
Additionally, three 843 MHz observations performed at the Mo- 
longolo Observatory Synthesis Telescope (MOST) and reported 
by Hannikainnen et al. (1998) are also simultaneous with these 
RXTE observations. 

This paper is structured as follows. We discuss the spectral 
analysis of archival Advanced Satellite for Cosmology and As- 
trophysics (ASCA) data in section 2. We look for evidence of 
Fe lines in the data and we characterize the soft (< 1 keV) X- 
ray data. In section 3 we present the RXTE data. We first dis- 
cuss the All Sky Monitor (ASM) data, and then we discuss the 
pointed observations. We perform spectral analysis much akin 
to that which we considered for Cyg X-l (Dove et al. 1998). 
Here, however, we consider Advection Dominated Accretion 
Flow (ADAF) models as well by using the models described 
by Di Matteo et al. (1998). We discuss the implications of the 
simultaneous radio data in section 4. In section 5 we discuss 
the implications of the X-ray observations for theoretical mod- 
els. We summarize our results in section 6. We present timing 
analysis of the RXTE data in a companion paper (Nowak et al. 
1999, henceforth paper II). 

2. ARCHIVAL ASCA OBSERVATIONS 

The ASCA archives contain four observations of the 
GX 339-4 region, one of which was sufficiently offset (> 3') 
so as not to detect the source. In appendix B, we describe the 
methods that we used to extract, filter, and analyze these ASCA 
observations. A log of the observations is given in in Table 1 . 
To the best of our knowledge, an analysis of these observations 
has not been published previously, except for a power spectrum 
for one of the observations (date not given, Dobrinskaya et al. 
1997). As we will discuss further below, the inferred 3-9 keV 
fluxes for the observations (Table 1) are lower by a factor of 
two to ten than the fluxes of the RXTE observations discussed 
in §3. 

We chose to fit the ASCA data with a phenomenological 
model consisting of a multicolor blackbody spectrum plus a 
broken power law, considered with and without a narrow Gaus- 
sian line at « 6.4 keV. These fits are similar to those performed 
for ASCA observations of the hard state of Cyg X-l (Ebisawa 
et al. 1996), which shows evidence for a weak and narrow Fe 
line with equivalent width % 40 e V, as well as for a soft excess 
well-modeled as a multicolor blackbody with peak temperature 
« 150eV. 

The fits with the phenomenological models yield rang- 
ing from 0.98 to 1.4. The brightest data set showed the greatest 
evidence for structure beyond this simple model. A sample fit 
is shown in Figure 1. Note that the neutral hydrogen column 
was fixed to 6 x 10 21 cm“ 2 . Allowed to freely vary, the neutral 
hydrogen column tended to float between 4 and 8 x 10 21 cm -2 , 
depending upon what combination of phenomenological mod- 
els was chosen, with minimal changes in the x 2 of the fits. As- 
sociated with these changes in best fit neutral hydrogen column 
were £ ±30% changes of the best fit peak temperature of the 
multicolor blackbody and even larger changes (factors of % 3) 
in the best fit normalization of the multicolor blackbody com- 
ponent. We should thus associate systematic error bars with 
these two parameters that are somewhat larger than the statisti- 


cal error b irs presented in Table 2. 

All fits improved with the addition of a narrow Gaussian line. 
In all fits we fixed the line width to 0. 1 keV (see Ebisawa et al. 
1996; whe always found o < 0.1 keV in fits to ASCA data of 
Cyg X-l), and for the lowest flux data set we also fix the line 
energy to 6.4 keV. For the lowest flux data set the A\ 2 = 5.5 for 
one additional parameter. By the F-test (Bevington & Robin- 
son 1992), this is an improvement to the fit for one additional 
parameter at the 98% confidence level. The other two data sets 
show even more significant improvements to the \ 2 - The best 
fit line equivalent widths ranged from as 30 to 60 eV. There is no 
compelling evidence for a strong flux dependence to the equiv- 
alent width of the line. 

The trarsition to the bright, soft state typically occurs at 3- 
9 keV lum nosities £ KT 9 ergs cm -2 s" 1 ; i.e., a factor of two to 
ten brighter than these ASCA observations. Thus these obser- 
vations offer useful tests of ADAF models, which are hypothe- 
sized to be most relevant to low-luminosity, hard state systems 
(Narayan 1996; Esin, McClintock & Narayan 1997). ADAF 
models predict a detectable correlation between the tempera- 
ture of the soft excess, the strength of the Fe line, and the source 
luminosity. They hypothesize that the luminosity decay of BH 
transients is due, in part, to an increase of the radius at which the 
accretion f ow transits from cold, geometrically thin, and radia- 
tively efficient to hot, geometrically thick, and advective (Esin, 
McClintock & Narayan 1997). In some models, the transition 
radius can grow to as large as 0(1O 4 Rq)> where R G = GM/c 2 . 
(Such larg«: transition radii are not a strict requirement of ADAF 
models; ir §3.2,4 we show that somewhat smaller transition 
radii, » 200-400 Rq, are preferred for ADAF models of the 
RXTE data.) As discussed by Esin, McClintock & Narayan 
(1997), one then expects the peak temperature of the soft ex- 
cess to decrease below 150 eV and the equivalent width of any 
Fe line to decrease to values less than « 30 eV. 

The best fit equivalent widths found for GX 339-4 are greater 
than can bn accommodated in ADAF models with a large tran- 
sition radius, and they are also slightly larger than predicted by 
the ‘sphen +disk* corona models described in §3.2.3 (see also 
Dove et al 1997; Dove et al. 1998). These latter models have 
a similar geometry to the ADAF models, and they often posit a 
coronal radius & 100 Rq. Likewise, we do not detect any large 
decreases in the best-fit disk temperatures with decreasing lu- 
minosity. Although it is dangerous to make a one-to-one cor- 
respondent t between a phenomenological fit component and a 
physical component, these best-fit values are suggestive of, but 
not definiti ve proof of, temperatures hotter than can be accom- 
modated in models where cold, soft X-ray emitting material 
exists at very large radii. 

3. RXTE OBSERVATIONS 
3.1. The Monitoring Campaign 

To stud) the long-term behavior of GX 339^, and to place 
our pointed observations within the context of the overall be- 
havior of tie source, we used data from the All Sky Monitor 
(ASM) on RXTE. The ASM provides lightcurves in three en- 
ergy band;, 1.3-3.0keV, 3.0-5.0keV, and 5.0-12.2keV, typ- 
ically comisting of several 90s measurements per day (see 
Levine et al. 1996; Remillard & Levine 1997; Lochner & 
Remillard 1997). In Figure 2 we present the ASM data of 
GX 339-4 up until Truncated Julian Date (TJD) « 1000 (1998 
July 6). We also indicate the dates of our RXTE observa- 
tions, as well as the dates of ATCA and MOST radio obser- 
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Table 1 : Log of the ASCA observations. All observations were taken in Bright 1-CCD mode. SISO: filtered SISO count rate. 


Obs. 

Date 

Integration time 
(ksec) 

SISO 
(cts s" 1 ) 

3-9 keV Flux 
(10~ 9 ergs cm -2 s" 1 ) 

1 

1994 August 24 

15 

3.6 

0.11 

2 

1994 September 12 

17 

6.4 

0.19 

3 

1995 September 8 

30 

17.5 

0.63 


Table 2: Parameters for a multicolor blackbody plus broken power law plus Gaussian line fits to ASCA data. 7i n : peak multicolor blackbody 
temperature. Adbb* multicolor blackbody normalization. Tj, IY broken power law photon indices. E t>: break energy. Ab P i: Power law normalization 
(photons keV 1 cm 2 s 1 at 1 keV). E\\ line energy. A\: Line normalization (photons cm~ 2 s -1 in the line). EW: line equivalent width. Uncertainties 
are at the 90% confidence level for one interesting parameter (Ax 2 = 2.71). 


Date 

7i„ 

Adbb 

r, 

E b 

r 2 

Abpi 

Ei 

Ai 

EW 

X 2 /DoF 

xid 


(keV) 

(xlO 4 ) 


(keV) 


(XlO* 2 ) 

(keV) 

(x 1CT 4 ) 

(eV) 



1994 Aug. 24 

O 1 4+0-01 
U< ^-0.02 

<7 < 7 + 0.9 
**- 1.6 

1 *70+003 
y 0 -0.03 

* -4-0.5 

1 £<7+0.04 

Lo2_0 0 3 

A 9+0.1 




1500/1439 

1.04 

1994 Aug. 24 

n ia+o-oi 

u* 1 4-0.02 

<7 < 7 + 0.7 
*•*- 1.6 

1 tq+0-03 
/0 -0.03 

<5 ^40.5 
* ) - J -0.6 

1 £4+0.03 
iO< +-0.03 

4 940. 1 
4 -*-0.l 

6.4 

0-6!^ 

34+25 

•^-19 

1494/1438 

1.04 

1994 Sept. 12 

0.15!™! 

<7 c + 0.6 
*•^- 1.0 

i q 1 +0.02 

1 * ol -0.0l 

<5 O+0.2 

J.8^0 2 

1 C £40.03 

7 ^+ 0.1 
’ * ^- 0. 1 




1603/1621 

0.99 

1994 Sept. 12 

o.i5iS:S8 

IggjM 

1 Q 1 40.01 
^Ol-O.Ol 

<3 O+0. 1 

i cQ + 0.01 

n < 7 + 0.0 
' -*- 0.1 

£ 7 £+0.08 

i-«2S3 

56!, 26 

1580/1619 

0.98 

1995 Sept. 08 

0 1 Q + 0-00 

U * iy -O.0l 


1 q<5+0.02 

3-7!S:l 

1 £1740.01 
1 -OU_o Q2 

25.6^3 




2597/1838 

1.41 

1995 Sept. 08 

O iq+0.00 
^-O.OO 


1 ’ V :> - 0.00 

3.7!S:i 

1 £7+0.01 
1 *^-0.01 

25.61? 

£ ci + 0.07 

&A1-0.06 

a 7 + 0.6 

- 3 --5-O.S 

40!^ 

2523/1836 

1.37 



Fig- I- — GX 339-4 ASCA observation of 1994 August 24. with energy bins rebinned by a factor of 25 for clarity. Model and associated residuals 
(data/model) are for the best fit multicolor blackbody plus broken power law without a Gaussian line component. For clarity, only the SISO (circles) 
and GIS2 (triangles) data are shown. 
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vations (Fender et al. 1997; Corbel et al. 1997; Hannikainnen 
et al. 1998). We discuss the long timescale variability of this 
lightcurve in paper II . 

Based upon model fits to the observations of Grebenev et al. 
(1991) (low and high state), and Miyamoto et al. (1991) (very 
high state), we expect the different states of GX 339-4 to have 
ASM count rates as indicated in Table 3. The ASCA and RXTE 
observations discussed here are most characteristic of weak to 
average luminosity hard states. A log of the RXTE pointed ob- 
servations and the simultaneous radio observations is given in 
Table 4. The transition to a higher flux level that occurs at TJD 
ft 800 has a characteristic flux of a high state, but is not as soft 
in the 2-10keV bands as the above cited high and very high 
states. This might be an example of what Mendez & van der 
Klis (1997) refer to as an ‘intermediate state’ between hard and 
soft. The variations observed in both the ASM lightcurve (prior 
to TJD ft 800) and the pointed RXTE observations discussed 
below represent more than a factor five variation in observed 
flux. Comparable variations have been observed in the radio, 
and furthermore the radio lightcurves show evidence of a cor- 
relation with both the ASM and Burst and Transient Survey 
Explorer (BATSE) lightcurves (Hannikainnen et al. 1998). 

3.2. PC A and HEXTE Observations 

In this section we present the results from our analysis of the 
data from the two pointed instruments on RXTE: the Propor- 
tional Counter Array (PCA), and the High Energy X-ray Tim- 
ing Experiment (HEXTE). See appendix B for a description of 
the instruments and of the details of the data extraction and pro- 
cessing. 

As we show in appendix B, there is a difference in the power- 
law slopes obtained from an analysis of spectra of the Crab with 
both instruments upon RXTE. In order to minimize the im- 
pact of this difference in the instrumental calibration onto the 
data analysis, we primarily analyze the data from both instru- 
ments individually and use the difference in model parameters 
between instruments as a gauge of the systematic uncertainties. 
We do perform some joint analysis of PCA and HEXTE data 
using various reflection models. In the following sections we 
discuss in detail the implications of the calibration uncertainty 
for our analysis. 

For our analysis of the RXTE broad band spectrum, we 
used several different spectral models to cover the range of 
parameterizations currently used in the literature to descripe 
the spectra of BHC. As in the ASCA analysis (§2), we fixed 
A/h = 6x 10 21 cm 2 . In the remainder of this section we de- 
scribe the various models in detail. We first used the purely 
phenomenological exponentially cutoff power law and broken 
power law models were for a broad characterization of the data. 
The results of this modeling are given in §3.2.1, the results of 
the fitting are presented in Tables 5 and 6. We then applied 
the three more physical models that are currently discussed 
in the literature: reflection of an power law off an (ionized) 
accretion disk (§3.2.2 and Table 7), ‘sphere and disk’ corona 
Comptonization models (§3.2.3 and Table 8), and ADAF mod- 
els (§3.2.4). The latter are only applied to the unfolded data. 

For the different models, we show the residuals for each 
of the best fits to the data from Observations 5 and 7 in Fig- 
ures 5, 3, and 4 (except for the ADAF models, where we present 
Observations 1 and 5). We chose to present these former two 
observations because not only are they at extremes in terms of 
observed luminosity (Observation 5 is the faintest, and Obser- 
vation 7 is the second brightest, Table 4), but also because they 


show characteristically different timing properties (paper II). 

3.2.1. Phenomenological Models 

Results from the purely phenomenological fits to the data, 
i.e., the broken power law and exponentially cutoff power law, 
are presented in Tables 5 and 6. We see that a broken power law 
plus a Gaussian models the PCA data extremely well. The ex- 
tremely low xitd (0.15-0.32) indicates that to some extent we 
may be firing systematic features in the PCA response. The 
same may be true for the best fit parameters of the Fe line. The 
line widths (o ft 0.6 keV) and equivalent widths (ft 1 lOeV) 
are larger than for the ASCA observations, with the exception 
of Observation 5 which has line parameters comparable to the 
ASCA observations. As Observation 5 has the lowest count 
rate, it is more dominated by statistical errors and less domi- 
nated by systematic errors than the other observations. Even 
ignoring the possible systematic effects, however, we see that 
any observed line is narrower and weaker than is commonly 
observed in AGN. 

The ft 3-10keV spectral power-law slope is close to the 
‘canonical value’ of T = 1 .7; however, the PCA shows evidence 
for a hardening of this spectral slope above ft lOkeV. HEXTE 
data alone also show the high energy spectrum to be harder 
than the 3-10keV spectrum (Fig. 4 and Table 6). Note that 
the difference between the PCA and HEXTE photon indices 
is greater han the discrepancy between the PCA and HEXTE 
best fit Crib photon indices (appendix B), and therefore it is 
unlikely tc be a systematic effect. 

Such a hard HEXTE spectrum is consistent with previous 
observation by the Oriented Scintillation Spectrometer Exper- 
iment (OSSE) on board the Compton Gamma-ray Observatory 
(CGRO) (Grabelsky et al. 1995; see also Zdziarski et al. 1998). 
Grabelsky et al. (1995) found a slightly harder photon index of 
T = 0.88 and an exponential cutoff of £cu t = 68 keV, somewhat 
lower than observed here. Note, however, that the OSSE ob- 
servations extended to ft 500 keV as opposed to the ~ 1 lOkeV 
of our HEXTE observation. Therefore, the HEXTE data for 
GX 339-4 do not strongly constrain the exponential rollover, 
and slightly harder power laws with lower exponential cutoffs 
are permit ed. 

3.2.2. Reflection Models 

A spectral hardening above ft 7keV is the expected sig- 
nature of reflection of a hard continuum off of cold material 
(Magdziars & Zdziarski 1995). Ueda, Ebisawa & Done (1994) 
applied ref ection models to Ginga data of GX 339^4 and found 
strong evidence of reflection, whereas Grabelsky et al. (1995) 
found no evidence of reflection in OSSE data. Zdziarski et al. 
(1998) jointly fit these simultaneously observed data sets and 
found that reflection models, albeit with a large Fe abundance, 
provide a very good description of the data. We have applied 
the modek of Magdziarz & Zdziarski (1995), as implemented 
in XSPEC (pexrav , pexriv ), to the GX 339-4 data. These mod- 
els consider an exponentially cutoff power law reflected off of 
neutral (pexrav ) or partially ionized ( pexriv ) cold material. 

In Table 7 we show the fit results for reflection off of partially 
ionized material similar to the models presented by Zdziarski 
et al. (199o)- Specifically, we fix the reflector inclination angle 
at 45°, fix he disk temperature at Ti* = 10 6 K, freeze the abun- 
dances at solar, but let the Fe abundance be a free parameter. In 
all our fits we found that the Gaussian line width, cr, would tend 
to drift toward 0, so we fixed o = 0.1 keV. For the combined 
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Fig. 2. — Left: RXTE All Sky Monitor data for GX 339-4 (5 day averages in the 1.3-1 2.2 keV band) vs. Truncated Julian Date (TJD) = Julian 
Date (JD) -2450000.5. Dashes indicate dates of our RXTE pointed observations, diamonds indicate dates of MOST radio observations (Han- 
nikainnen et al. 1998), and triangles indicate dates of ATCA radio observations (Fender et al. 1997, Corbel et al. 1998). Right: ASM colors 
for GX 339-4 vs TJD. Colors shown are: (1.3-3 keV)/(3-5 keV), (3-5 keV)/(5-12 keV)+5, and (1.3-3 keV)/(5-12.2 keV)+10 (the latter two colors 
have been offset for clarity). 


Table 3: Approximate expected ASM colors for the different states of GX 339—4. See text for explanation. 


State 

count rate 
ctss” 1 

1.3-3 keV/ 
3-5 keV 

1.3-3keV/ 

5-12.2keV 

low 

7 

2 

2 

high 

15 

4 

30 

very high 

60 

4 

15 


Table 4: Log of RXTE and radio observations. We list: T, the duration of the RXTE observations; the average PC A count rate; the average 
(3-9 keV), (9-30 keV), and (30-100 keV) energy fluxes (all normalized to the PCA calibration); the flux density of the 8.3-9. 1 GHz ATCA obser- 
vations; the flux density of the 843 MHz MOST observations; and a = A logy/ A logS Vi the spectral index of the radio observations (Fender et al. 
1997; Corbel et al. 1997; Hannikainnen et al. 1998). 


Obs 

Date 

TJD 

T 

(ksec) 

Rate 

(cps) 

3-9 keV 

9-30 keV 
(10" 9 erg cm -2 s -1 ) 

30-100 keV 

ATCA 

(mJy) 

MOST 

(mJy) 

Q 

01 

1997 Feb. 02 

481 

11 

830 

1.07 

1.68 

2.65 

9.1 ± 0.2 

7.0 ±0.7 

0.11 

02 

1997 Feb. 10 

489 

10 

730 

0.94 

1.50 

2.41 

8.2± 0.2 

6.3±0.7 

0.11 

03 

1997 Feb. 17 

496 

8 

700 

0.90 

1.43 

2.35 

8.7 ±0.2 

6.1 ±0.7 

0.15 

04 

1997 Apr. 29 

567 

10 

470 

0.60 

0.97 

1.55 




05 

1997 Jul. 07 

636 

10 

200 

0.25 

0.43 

0.75 




06 

1997 Aug. 23 

683 

11 

650 

0.74 

1.18 

1.98 




07 

1997 Sep. 19 

710 

10 

730 

0.96 

1.48 

2.36 




08 

1997 Oct. 28 

749 

10 

480 

0.63 

1.01 

1.68 
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Fig. 3. — Spectral modeling of the PCA data from Observation 5 (left) and 7 (right). Residues are shown as the contribution to x- a ) Count rate 
spectrum and the best-fit broken power-law with Gaussian line, b) Contribution to x from the broken power-law with Gaussian line, c) Contribution 
to x from the ionized reflection model (pexriv) with Gaussian line, d) Contribution to x f om the sphere+disk model with additional Gaussian line. 


Table 5: Parameters for Gaussian line plus broken power law models (PCA only). 


Model: const. * highecut^cu^Efow) exp(-/V , Ha fl )(gauss(£u 1Kj cr,Ab*)+bknpower(ri ) r 2 1 £bf«k,Abkn)) 


ObS £\mc G Aline EW T 1 £ break T2 Abkn X J /d0f xL 

(keV) (keV) (keV) (keV) (xl0‘ 3 ) (eV) (keV) 


01 

6.483!;!; 

0.631 

1.973£ 

1 30*13 

i.8o3:g! 

11.235 

1 C'j+O-M 

0.443°! 

13.6/ 52 

0.26 

02 

6.47t!’ 

0.631 

1.7333! 

12935° 

180383 

10.93; 


0.383;“! 

16.2/52 

0.31 

03 

6.4m; 

o.53:l 

i.553:ll 

12131° 

1-793:8! 

10.93J 

1 ct+O-02 

1 * J J-0,03 

o.373;8! 

12.8/52 

0.25 

04 

6.483;“ 

o.53:l 

o.993;2 

11635 s 

1773:?! 

10.83:5 

1 

1 * JJ -0.03 

0.233?? 

14.7/52 

0.28 

05 

6.473;!| 

0.231 

0.273:® 

7433 

i.7o3?! 

10.83! 

i-483;2 

0.083:2? 

17.7/52 

0.34 

06 

6.433!! 

o.63:l 

1183® 

11331 6 

1-773:?! 

10.93:5 

1-503:® 

n oq+O.oo 
u - zo -0.00 

13.2/52 

0.25 

07 

6.503!! 

0.631 

1.743 6 3« 

1313 s0 

1.803?! 

10.935 

1 c-i+0‘02 

0.393;°! 

12.8/52 

0.25 

08 

6.453:!! 

o.63:l 

1.083:1! 

12135° 

1-763?! 

10.83 s 

i r a +0.03 
1 *- )J -0.03 

0.24388 

15.2/52 

0.29 
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Fig. 4. — Spectral modeling of the HEXTE data from observation 5 (left) and 7 (right). Residues are shown as the contribution to x ■ a) Count rate 
spectrum and the best-fit power-law with exponential cutoff, b) Contibution to \ from the power-law with exponential cutoff, c) Contribution to \ 
from the sphere+disk model. 


Table 6: Parameters for exponentially cutoff power law models (HEXTE only). 

Model: const. (i4 P L£ r ’ r exp(--£/£ CU i)) 


Obs const. T £ cut Apl const. x 2// dof xLj 

(keV) 


01 

LOO 

1 o«+0-06 

101*12 

0.082tm! 

n QQ+ 0-01 

69.2/ 80 

0.87 

02 

LOO 

1 t 'Y*Q08 
1 • 1 *-0.08 

79112 

0-052^ 

0.92t^ 

69.0/ 80 

0.86 

03 

LOO 

1 16-0.18 

94 12 . 

0.052!™?, 

1.08t“J 

71.8/ 80 

0.90 

04 

LOO 



0.035^ 

0 98!^;* 

76.2/ 80 

0.95 

05 

LOO 

1 IR 4 * 319 
• • t ®-0.23 

115!“ 

0.016^ 

0.99t« 

71.3/ 80 

0.89 

06 

LOO 

i i q+0 09 

103!l? 

0.049!S;81g 

0.93!{im 

87.9/ 80 

1.10 

07 

LOO 

1 0 1 +0-07 

i-^-o.os 

95!® 

0 . 066 !™“ 

0.961S 

101.5/ 80 

1.27 

08 

LOO 

b 

oo 

oo 

81!!? 

0.03l!™°? 

0.98!&{g 

106.6/ 80 

1.33 








normalized counts/sec/keV 
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Fig. 5. — Spectral modeling of the joint PCA and HEX I t data for Observation 5 and 7. Residues are shown as the contribution to x • a ) Count 
rate spectrum and the best fit broken power-law with Gaussian line, b) Contribution to y from the best fit broken power-law with a Gaussian line, 
c) Contribution to x from the best fit ionized reflector model with a Gaussian line. 
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Table 7: Parameters for Gaussian line plus ionized reflection models after Magdziarz & Zdziarski (1995). Models have been fit to: PCA data only 
(52 dof); PCA+HEXTE data where the photon index of the reflected power law has been constrained to be identical for the PCA and the HEXTE 
data (134 dof); and PCA+HEXTE data where the photon index of the reflected power law has been allowed to differ between the PCA and the 
HEXTE data (133 dof). Fit parameters are the power law slope, T; cutoff energy, £f 0 i<j; relative reflection fraction, / = Afy27r; Fe abundance 
relative to solar, A?c\ the disk ionization parameter, £ = luminosity/(density x radius 2 ); and a normalization, A^ x . The ionized reflection models 
also have parameters for the temperature of the disk, 7duk (fixed at 10 6 K); abundances of elements heavier than He relative to solar, Ax (fixed at 
unity); and disk inclination angle, i (fixed at 45°). Parameters typeset in italics have been held constant for that particular fit. The iron line width 
was held fixed at 0. 1 ke V. 
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PCA and HEXTE data, we also fixed the Gaussian line energy 
to 6.4 keV. For fits to PCA data alone and joint PCA/HEXTE 
data, the exponential cutoff energy, £f 0 id» would drift towards 
very large energy (^> 1000 keV). We therefore only considered 
pure power laws without cutoffs. Zdziarski et al. (1998) have 
argued that the high energy cutoff is sharper than exponential, 
which one would not expect to be strongly constrained by the 
combined PCA/HEXTE data. 

As for the Ginga data of GX 339-4 (Ueda, Ebisawa & Done 
1994), the PCA data alone were extremely well described by 
reflection models. Again, however, the extremely low x^ (as 
low as 0.16) makes us caution that these models might partly 
be fitting systematics in the PCA response. PCA and Ginga 
are also very similar instruments in terms of design, and so to 
some extent they should exhibit similar systematic effects (As 
discussed in appendix B, the internal consistency of the PCA 
calibration is now as good as or better than that for the Ginga 
calibration.). Note that the best fit Fe line equivalent widths 
here are significantly smaller than those found with the purely 
phenomenological models discussed in §3.2.1 . 

The fits for the HEXTE data alone (not shown) were simi- 
lar to the OSSE results of Grabelsky et al. (1995). Namely, if 
one allows for an exponential cutoff (typically « lOOkeV) to 
the power law, the best fit reflection fraction becomes / £ 0.01 . 
Such a small reflection fraction is not surprising considering 
how well a pure exponentially cutoff power law fits the HEXTE 
data (Table 6). If one does not allow an exponential cutoff, the 
reflection fraction becomes / £ 3. Such a fit is trying to mimic 
a hard power law with a high energy cutoff. 

A joint analysis of the PCA and the HEXTE data should be 
similar to a joint analysis of the Ginga and OSSE data. Indeed, 
such an analysis yields results comparable to those presented 
by Zdziarski et al. (1998) if we constrain the photon index of 
the incident power law to be the same for both the PCA and 
HEXTE data. Notable in these results is the fairly large over- 
abundance of Fe (4p c « 4). Similarly, Zdziarski et al. (1998) 
find a large Ap e = 2. 5-3.0 except for a short data set, more 
likely dominated by statistical errors rather than systematic er- 
rors, where they find Af* = 1 .6-2.0. For our joint PCA/HEXTE 
data, the best-fit reflection fraction is approximately 20% larger 
than the best-fit reflection fraction for PCA data alone. Such 
an increase in reflection fraction in general will reproduce the 
spectral hardening seen in the HEXTE energy bands. Increas- 
ing the average best fit Fe abundance from (Ape) = 14 (PCA 
data alone) to (Ap*) = 4. 1 (joint PCA/HEXTE data) also leads 
to an increased spectral hardening above « 7 keV, while leav- 
ing the spectrum below « 7 keV relatively unchanged. That is, 
such a fit helps to reproduce the spectral break at « 10 keV. 

For the joint PCA/HEXTE analysis, there is clearly a worry 
that these results arc influenced by the systematic differences 
between the PCA and the HEXTE responses. We therefore per- 
formed reflection model fits where we constrained all fit param- 
eters to be the same for the PCA and the HEXTE data except for 
the incident power law photon index, which we allowed to vary 
between the two instruments 1 * . Such models provided reason- 
ably good fits to the data, with x„d ranging from 0.72 to 1.25. 
The difference between the PCA and the HEXTE best fit pho- 
ton indices ranged from 0.08 to 0.16, with an average value of 
0.12. This is consistent with the systematic difference between 
the best-fit photon indeces for the Crab spectrum. For these 


models we find (Ap*) = 1.9, which is more consistent with the 
results for fits to the PCA data only, and is slightly smaller than 
the results found by Zdziarski et al. (1998). Note that we also 
find smaller values of the ionization parameter, £, than were 
found by Zdziarski et al. (1998). 

3.2.3. Corona Models 

We considered 4 sphere+disk’ Comptonization models (Dove 
et al. 1997) of the GX 339-4 observations. We have previously 
applied these models to an RXTE observation of Cygnus X-l 
(Dove et al. 1998). The models consist of a central, spheri- 
cal corona surrounded by a geometrically thin, flat disk. Seed 
photons fcr Comptonization come from the disk, which has a 
radial temperature distribution £7di S k(#) oc R ~ y4 and a temper- 
ature of 150eV at the inner edge of the disk. Hard flux from 
the corona further leads to reflection features from the disk or 
to soft photons due to thcrmalization of the hard radiation. The 
(non-uniform) temperature and pair balance within the corona 
are self-consistently calculated from the radiation field (Dove 
et al. 1997). 

As described by Dove et al. (1997), we parameterize our 
models by the coronal compactness 


L = aT 


m^c 3 R c 


( 1 ) 


where cj is the Thomson cross section, is the electron 
mass, Lq is the luminosity of the corona, and Rc is the ra- 
dius of the corona. Likewise, we can define a disk compact- 
ness, £d = (1 -fc)( 0 VmtC*)Pc/Rc> where Pq is the total rate 
of gravitational energy dissipated in the system, and f c is the 
fraction dissipated in the corona. In calculating the numerical 
models, we set l& = 1 . Models with other values of £4 yield the 
same ranges of self-consistent coronal temperatures and opac- 
ities, and in general f c =£/(£ d+4) (Dove et al. 1997). Based 
upon the 4 sphere+disk’ geometry, a fraction / « 0.32 of the 
coronal flux is absorbed by the disk (Dove et al. 1997). The 
models are further parameterized by an initial electron coronal 
optical depth, r c (approximately equal to the total optical depth, 
as pair production is negligible for the parameters of interest to 
us), and a normalization constant Ak 0{ . From the best-fit com- 
pactness and optical depth, the average coronal temperature can 
be calculated a posteriori. 

Attempts to fit these models to the joint PCA/HEXTE data 
failed. Topical xLi values, even allowing for the inclusion 
of an extra Gaussian line component, were £ 1.3. These fits 
showed a dear tendency for a hardening in the HEXTE band, 
and therefi ire we considered them to be influenced by the cross- 
calibration uncertainties between the PCA and the HEXTE in- 
struments (Note that our previous fits of Cyg X-l used an 
earlier version of the PCA response where we applied 1.5% 
systematic uncertainties across the entire PCA band; these fits 
yielded Xnd ^ 1*6, without considering an additional Gaussian 
componen:.). We therefore considered ‘sphere+disk’ models fit 
to the PCA and the HEXTE data separately. In Table 8, we 
present the best-fit parameters for these models applied to our 
GX 339-4 data. 

Although our numerical ‘sphere+disk’ models do include 
reflection uid a flourescent Fe line (typical equivalent width 
« 25 cV) from the disk, the PCA data showed residuals in the 


1 The photon index was constrained to be the same for HEXTE Ouster A and B. The necessary different normalizations between the PCA and the HEXTE models 

were subsumed into the constants multiplying the HEXTE models. As HEXTE requires a hardet power law, these constants were now 0.39-0.52, as opposed to 

~ 0.67. Furthermore, the constants now show larger uncertainties, as the uncertainty of the HEXTE photon index now couples strongly to the value of the constants. 
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Table 8: Models of ‘sphere+disk’ Comptonization plus a Gaussian line fit to PCA data only. The three fit parameters of the Comptonization model 
are the compactness of the corona, £ c , the coronal optical depth, r c , and a normalization constant, A kot- The Gaussian line is parameterized as in the 
previous tables. From the best fit parameters, the equivalent width of the line, EW, and the density averaged coronal temperature, kT c , are derived. 


Obs 

feline 

O 

Aune 

EW 

1 c 

T c 

Aim 

const. 


X 2 /dof 

xL 


(keV) 

(keV) 

(xl0~ 3 ) 

(eV) 





(keV) 



01 

6.39!S:!? 

0.83? 

7 OQ + 0 - 60 
J - w -0.53 

197!?“ 

r\ £7+0.05 
V.O 2 _q 0 4 

-a 7Q+O.06 

7 70 +O.O 6 
°-0.08 


28.631 

19.2/53 

0.36 

01 


... 



1 Q7+O.I7 
l* w - 0.13 

7 0+0.5 
''-0.3 

A Q 1 +0-09 
UV1 -0.09 

0.9933 

43.9$ 

71.2/80 

0.89 

02 

6-36^;l| 

0.83;? 

7 AQ 40 - 57 
zoy -0.50 

194!8 

0. 03^)05 

•a 77+0.07 

J. JJ _0 07 

7 rj 7 +O.O 6 
ZU/ -0.07 


27.93? 

15.7/53 

0.30 

02 





1 QA+012 
I j 2 

7 A+O - 3 
3*0^0 3 

A 77+0.04 
U /: >-0.03 

0.9233 

34.0!? 

67.6/80 

0.84 

03 

6.363$ 

o.83j 

9 47+0+4 

182!?? 

£/:+0.06 

U.OO- 0.06 

7 7 <+0.06 
3*3j_oo7 

1.953$ 


27.931 

16.7/53 

0.31 

03 





O 0 1 +0 *48 

3.431 

0.613$ 

10833 

37.4!? 6 

72.4/80 

0.90 

04 

6-40!£n 

o.83:? 

i.533:g 

174!?, 

0 . 703 $ 

-a cc+0.08 

J. JJ_0,09 

1.273$ 


26.631 

13.5/53 

0.25 

04 


... 



2093?? 

3.43? 

0463$ 

0.9833 

36.93 

78.1/80 

0.98 

05 

f. 70 +O .31 

O.ZO_o97 

0.93J 

0.523:?? 

1 70+100 
1 J O-I00 

fl*h 8 ^oB 

4.393?? 

0.5333 


20.631 

20.5/53 

0.39 

05 


• •• . 



^•*> y -0.65 

3-l!?i 

A 1 a+O.03 

U* 1 0-0.03 

0.9933 

44.9!l? 3 

71.8/80 

0.90 

06 

6.263:8 

0.93:? 

7 AA+0.50 

1863? 

A AA+O-07 
O.00_ox)5 

3603$ 

1 £7+0.04 
1 •*> '- 0.06 


26.231 

16.5/53 

0.31 

06 


... 



7 77 + 0.21 
z - j:> -0.18 

7 1 + 0.6 
^ * * -0.4 

A c^+0.08 
U:>4 -0.04 

0.9333 

42.3!? 1 

88.7/80 

1.11 

07 

6.403:1? 

0.83:? 

7 cs 7+0.56 
Z.D7_o49 

192!?? 

A A^+O-06 

O.Oj«oo5 

7 70 +O.O 7 
■*•■**>-0.06 

7 rv/;+0.05 
2 . 00 _oq 7 


28.431 

15.6/53 

0.30 

07 





1 QQ+013 
1<w -0.16 

7 1+0.4 
• 31 -0.4 

A 7«+0.00 
7 -’-0.04 

0.9633 

40.73° 

102.9/80 

1.29 

08 

6.3438 

o.93:? 

1 74 + 0+0 
t* '^-0.40 

1893? 

o.7i3$ 

7 ^7+0.09 
3*t>/«o,09 

1 70+0-05 


25.931 

15.5/53 

0.29 

08 





2.373:?? 

3731 

0 4 7+0.03 

0.9833 

34.6!? 

109.0/80 

1.36 


5-7 keV band, similar as in in our fits to RXTE data of Cyg X-l 
(Dove et al. 1998). We included an additional Gaussian com- 
ponent to our fits. The equivalent widths of the additional lines 
were « 150eV, and they appeared to be broad (a & 0.8 keV). 
This additional line may be attributable partly to remaining un- 
certainties in the PCA response. For these fits, as well as for 
the reflection model fits, lines with energies significantly less 
than 6.4 keV can be fit, and this is likely a systematic effect. 
Part of the discrepancy between the data and the model, how- 
ever, is significant. As we have discussed for our fits to the 
RXTE Cyg X-l data, there are several possible physical inter- 
pretations for the additional required equivalent width: There 
may be an overlap between the disk and the sphere (our models 
invoke a sharp transition), the disk may be flared (we model 
a flat disk), the disk may have non-solar abundances, or al- 
ternatively one might invoke a ‘patchy disk* embedded in the 
corona (Zdziarski et al. 1998). The best-fit reflection fractions 
of / « 0.4-0.5 found above are further indication that our mod- 
els may require an additional source of reflected flux. 

Allowing an additional Gaussian line component, the fits to 
the PCA data yield extremely low which could be indi- 
cating that we are partly fitting systematic features in the PCA 
response. Note also that the PCA data fits yielded consistently 
larger optical depths and consistently lower compactness pa- 
rameters than the HEXTE data fits. The latter was more sig- 
nificant, and is again indicative of the HEXTE response being 
harder than the PCA response. Both instruments yielded op- 
tical depths r c « 3-4; however, due to the discrepancy in the 
best-fit spectral slopes between the PCA and the HEXTE bands, 
the best-fit average coronal temperatures range from 21-30keV 
(PCA) to 34-45 keV (HEXTE). 

3.2.4. ADAF Models 

The basic picture of mass accretion via an ADAF in the con- 
text of galactic BHC was introduced by Ichimaru (1977) and 


has been elaborated upon in a series of papers by Narayan and 
collaborators (Narayan, Kato & Honma 1997; Esin, McClin- 
tock & Narayan 1997). The accretion flow is divided into two 
distinct zones: the inner part is modeled as a hot, optically thin 
ADAF similar in some respects to the spherical corona dis- 
cussed above, while the outer part consists of a standard op- 
tically thick, geometrically thin disk. The transition radius be- 
tween the two zones, r* = R t /Ro y is one of the model parame- 
ters. We compute the ADAF spectrum according to the proce- 
dure described by Di Matteo et al. (1998). The electrons in an 
ADAF cool via three processes: bremsstrahlung, synchrotron 
radiation and inverse Compton scattering. In addition we add 
the emission from a thin disk — calculated as a standard mul- 
ticolor blackbody — and include the Compton reflection com- 
ponent due to the scattering of high energy photons incident on 
the disk. 

In the ADAF models discussed here, we fix the black hole 
mass to be m = M/M 0 = 6, assume the magnetic field to be in 
equipartition with thermal pressure ((3 - 0.5), and set the stan- 
dard Shakura-Sunyaev viscosity parameter (Shakura & Sun- 
yaev 1973) to be ass = 0.3. We normalize the accretion rate 
to m = McVLedd, where L Edd is the Eddington luminosity of 
the source. The hard state corresponds to mass accretion rates 
m < merit = 10“ 2 , where rhan is the critical value above which 
an ADAF no longer exists. As m increases towards m cr i t , the 
scattering optical depth of the ADAF goes up which causes the 
spectrum to become harder and smoother. Most of the flux from 
the ADAF plus disk configuration is emitted around lOOkeV 
and the spectrum falls off at higher energies. 

The model spectrum changes mainly as a function of r 0 and 
m. The various spectral states correspond to different values of 
these parameters. For example, Esin, McClintock & Narayan 
(1997) attempt to explain the initial transition from soft to hard 
seen in the decay of Nova Muscae by a large change in r n (from 
rtr « 10 to rtr « 10 4 ), followed by an exponential decay in m for 
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the subsequent evolution of this transient system. The ASCA 
data of GX 339-4 discussed in §2 imply that comparably large 
changes in r n are not relevant to those observations. Here, how- 
ever, unfolded RXTE data from Observation 1 and Observa- 
tion 5, the brightest and faintest observation respectively, can 
be described by ADAF models with r * = 200, m = 0.08 and 
rtr = 400, m = 0.05, respectively. These model spectra and 
RXTE spectra for Observations 1 and 5 unfolded with a cut- 
off broken power law plus Gaussian line are shown in Fig. 6. 
In these ADAF models, the observed spectral and luminosity 
changes of GX 339-4 are predominantly driven by changes in 
the transition radius; the implied accretion rate change is sub- 
stantially smaller than the factor of 5 for the observed luminos- 
ity change. 

For r u ~ 200, which provides a rough description of Obser- 
vation 1, the disk blackbody emission peaks in the far UV/soft 
X-rays and dominates over the synchrotron emission, which 
peaks in the optical/UV (see also Zdziarski et al. 1998). The 
peak synchrotron emission frequency scales as oc m~ v2 m u2 i% 4 
and peaks in the range v = 10 n -10 12 Hz for supermassive black 
holes and v = 10 15 -10 16 Hz for galactic black holes. The spec- 
trum below the peak is approximately S^oci/ 2 . The synchrotron 
emission can contribute significantly to the radio emission of 
super-massive black holes (although see Di Matteo et al. 1998); 
however, the predicted radio flux of GX 339—4 is ten orders of 
magnitude below the observed 7 mJy flux at 843 MHz. Thus, 
there must be an extended source of radio emission, which we 
further discuss in the next section. 


4. SIMULTANEOUS RADIO OBSERVATIONS 

The first three of our RXTE observations were simultane- 
ous with 843 MHz observations taken with the Molongolo Ob- 
servatory Synthesis Telescope (MOST), and with 8.3-9. 1 GHz 
observations taken at the Australian Telescope Compact Array 
(ATCA). Extensive discussion of the MOST and ATCA obser- 
vations can be found in Hannikainncn et al. ( 1 998), Corbel et al. 
(1997), respectively (see also Table 4). 

An estimate of the minimum size of the radio emitting re- 
gion can be obtained by noting that observationally the bright- 
ness temperatures of radio sources usually are not larger than 
10 12 K, else the electrons will suffer catastrophic inverse Comp- 
ton losses. The brightness temperature of a uniformly bright 
spherical source is given by ( cD/dv^S^/lnk , where d is the di- 
ameter of the source, D is its distance, u is the observed radio 
frequency, S v is the observed radio flux density, c is the speed 
of light, and k is the Boltzmann constant. Taking the 7 mJy 
observed at 843 Mz by MOST and the fact that GX 339-4 is 
unresolved, we derive 
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which is orders of magnitude larger than the inferred size of 
the X-ray emitting region, even for models that posit extremely 
extended coronae (e.g., Esin, McClintock & Narayan 1997; 
Kazanas, Hua & Titarchuk 1997). 

This size scale is oc i/" 1 , so that emission at 8.6 GHz could 
arise in a region an order of magnitude smaller than that respon- 
sible for the emission at 843 MHz. Indeed, the flat spectrum 
emission is likely to arise in a conical jet with a radially de- 
creasing optical depth (e.g., Hjellming & Johnston 1988). Thus 
the outflow likely has an extent of O(10 7 GM/c 2 ) or greater. 
Similar estimates for source size have been made for the other 


persistent black hole candidate and Z-source neutron star X-ray 
binaries by Fender & Hendry (1999). 

Assuming a radio spectral index of a = 0.1 and a sharp cut- 
off at 10/iin (a reasonable upper limit for where the radio flux 
becomes optically thin, and typical of where ADAF models be- 
come optically thin in the radio), the radio flux is approximately 
0.1% of the 3-100keV X-ray flux. The correlation between the 
X-ray and radio fluxes found by Hannikainnen et al. (1998), 
comparable to the X-ray/radio correlation observed in Cyg X- 
1 (Pooley, Fender & Brocksopp 1998), suggests that there is a 
coupling between the inner accretion disk and the extended out- 
flow on timescales of 7 days or less. Matter leaving the corona 
at the escape velocity (0.25c at 30 Rq) and thereafter decelerat- 
ing under the influence of gravity would take roughly 7 days to 
travel a distance of 10 7 Rq. As 7 days is the upper bound to the 
correlation timescale, the radio emitting outflow must leave at 
slightly greater than escape velocity, or there must be at least 
some amount of acceleration of the outflow. 

Although the radio observations are strictly simultaneous 
with our first three RXTE observations, GX 339-4 exhibits less 
than 1% rms variability over the shortest timescales for which 
a reasonable radio flux estimate can be made (> 10 minutes). 
Thus there are no strong features to correlate between the radio 
and X-ray bands. 


5. DISCUSSION 

Coronal She and Luminosity Variation : 

The relationship between the inferred size of the corona and 
the magnitude of the observed flux depends upon which spec- 
tral model we are considering. As discussed in §3.2.4, for 
ADAF models one can associate lower fluxes with increased 
coronal radii. A larger coronal radius implies a lower efficiency 
and hence a decreased observed flux, even for constant accre- 
tion rates. Paper II shows that the characteristic power spectral 
density (PSD) timescale for GX 339—4 decreases for the low- 
est observed flux (Observation 5). If one associates the PSD 
timescale with characteristic disk timescales, this could be in 
agreement with an increased coronal radius. However, in pa- 
per II we also show that the time lags between hard and soft 
X-ray variability decreases with decreasing flux, which seems 
counter to i positive correlation between flux and coronal size. 

The ‘spiere+disk* coronal models make no assumptions 
about the radiative efficiency of the accretion. The flux can 
be either positively or negatively correlated with coronal ra- 
dius, deper ding upon the variations of the coronal compactness, 
4, and the temperature, 7i, at the inner edge of the accretion 
disk that s irrounds the corona (Dove et al. 1997). Note that 
the ‘sphere+disk’ models used in this work, contrary to many 
ADAF models, do not consider synchrotron photons as a source 
of seed photons for Comptonization. 

Using the definitions of £ c , ^ d , and / given in §3.2.3, energy 
balance in Jie ‘sphere+disk* system determines the coronal ra- 
dius, to within factors of order unity, to be given by 



where M is the mass of the compact object, D is the distance to 
the source, and F w is the bolometric flux of the source. If r d , 
/, and t c were held fixed, then the coronal radius would be 
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Fig. 6. — Advection Dominated Accretion Flow models for the unfolded RXTE data from Observation 1 (solid line) and Observation 5 (dash-dot 
line). Parameters consistent with the unfolded data are described in the text. A source distance of 4 kpc was assumed. 


positively correlated with flux. Whereas this might pose some 
problems for understanding the flux dependence of the charac- 
teristic timescales observed in the PSD, this would agree with 
the flux dependence of the X-ray variability time lags (paper II). 
However, as the RXTE bandpass does not usefully extend be- 
low « 3 keV, we do not have a good handle on the flux depen- 
dence of Ta. If 7d a F t f t with (3 > 1/4, then increasing flux 
could imply decreasing coronal radius. ASCA observations si- 
multaneous with a campaign of RXTE observations could help 
answer this question. 

Correlations Among Spectral Parameters: 

Ueda, Ebisawa & Done (1994) claimed that reflection mod- 
els of GX 339-4 exhibited a correlation between photon in- 
dex, T, and reflection fraction, /, with softer spectra implying 
greater reflection. Zdziarski (1998) has claimed that this corre- 
lation extends to reflection models of Seyfert 1 galaxies as well. 
Such a correlation is not unreasonable to expect. For example, 
if we allow the corona and disk to overlap to some extent in 
the ‘sphere+disk’ model (Poutanen, Krolik & Ryde 1997), then 
we expect the increase in the flux of seed photons to cool the 
corona and lead to a softer spectral index. Likewise, the cover- 
ing fraction of the disk would be increased, in agreement with 
the suggested correlation. 

In Figure 7 we plot / vs. T for our reflection model fits to 
GX 339-4. Contrary to the claims of Ueda, Ebisawa & Done 
(1994) and Zdziarski (1998), however, there is no strong ev- 
idence for a correlation. Fitting the reflection fraction with a 
function linear in T, as opposed to fitting with the mean value 
of /, improves the x 2 of the fits by 0.2 to 0.8, which is not 
significant. Fitting with the mean gives X'Li = 0*2 t0 1 -8. 

We do note two possible trends from the reflection model 
fits. First, as has been noted for other hard state galactic black 
hole candidates (Tanaka & Lewin 1995 and references therein), 
there may be a correlation between flux and photon index with 
lower flux implying a harder source. (The significance of the 
correlation is driven by Observation 5, the faintest and hard- 


est of the observations. However, a similar correlation is also 
present in color-intensity diagrams.) Such a correlation is con- 
sistent with the expectations of ADAF models where the ra- 
dius increases with decreasing flux (Figure 6). Again, our 
‘sphere+disk* corona models do not predict a clear trend with- 
out knowing the flux dependences of other parameters such as 
T*. 

Second, the ionization parameter, £, is positively correlated 
with flux. Such a correlation was noted by Zycki, Done & 
Smith (1998) for Ginga observations of Nova Muscae. It is 
not unreasonable to expect the disk to become increasingly ion- 
ized with increasing flux. We again caution that it is danger- 
ous to make one-to-one correlations between a model fit pa- 
rameter and a true physical parameter. Furthermore, the sig- 
nificance of the correlation is again almost entirely determined 
by Observation 5, the faintest observation, which has £ « 0. 
However, if we take the flux dependence of £ as being real 
and interpret it physically, it provides some constraints on the 
flux dependence of the coronal radius. The ionization param- 
eter is oc Fto/(pR 2 ), where p is the density of the disk. For 
a gas pressure-dominated Shakura-Sunyaev a-disk, p a /C 1 65 
(Shakura & Sunyaev 1973). In order for £ to be roughly lin- 
ear in flux (the actual dependence is not strongly constrained 
by the data), we require Fu* F 1 35 . (For our ‘sphere+disk* 
models this would further require g$ F" 1 ' 3 , depending upon 
the flux dependences of £ c , g, etc.) Taken physically and 

in the context of a gas pressure dominated Shakura-Sunyaev a- 
model, the flux dependence of £ implies that the coronal radius 
increases with increasing flux. 

6. SUMMARY 

We have presented a series of observations of the black hole 
candidate GX 339—4 in low luminosity, spectrally hard states. 
These observations consisted of three separate archival ASCA 
and eight separate RXTE data sets. All of these observations 
exhibited (3-9 keV) flux < 10 -9 ergs s“ ! cm -2 , and the observed 
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Fig. 7. — Left: Reflection /racti on vs. photon index, I\ for models fit to PCA data only (squares), models fitting the same photon index to PCA 
and HEXTE data simultaneously (triangles), and models allowing the PCA and HEXTE j hoton indeces and normalizations to be different (squares; 
HEXTE photon index shown). Middle: Photon index, I\ vs. observed 3-9 keV flux for the same reflection models as on the left. Right: Disk 
ionization parameter, f , in units of ergs cm s -1 for the same reflection models as on the 1< ft. 


fluxes spanned roughly a factor of 5 in range for both the ASCA 
and RXTE data sets. Subject to uncertainties in the cross cali- 
bration between ASCA and RXTE, the faintest ASCA observa- 
tion was approximately a factor of two fainter than the faintest 
RXTE observation. 

All of these observations showed evidence for an « 6.4 keV 
Fe line with equivalent widths in the range of « 20-1 00 eV. 
The ASCA observations further showed evidence for a soft 
excess that was well-modeled by a power law plus a multi- 
color blackbody spectrum with peak temperatures in the range 
ss 150-200eV. Both of these factors considered together argue 
against ‘sphere-fdisk’ or ADAF type-geometry coronae with 
extremely large coronal radii of C?(10 4 Rq) (e.g., Esin, McClin- 
tock & Narayan 1997). 

The RXTE data sets were well-fit by 4 sphere+disk’ Comp- 
tonization models with coronal temperatures in the range 20- 
50keV and optical depths in the range of r w 3. These fits were 
similar to our previous fits to RXTE data of Cyg X-l. Advec- 
tion Dominated Accretion Flow models, which posit a similar 
geometry, also provided reasonable descriptions of the unfolded 
RXTE data. These latter models can imply that the coronal ra- 
dius increases with decreasing flux. The 4 sphere+disk’ corona 
models do not make a specific prediction for the dependence 
of the coronal radius on the flux; however, they can be consis- 
tent with a positive correlation between coronal radius and flux. 
As described in paper 13, a postive correlation between flux and 
coronal radius is consistent with the observed flux dependence 
of the time lags between hard and soft X-ray variability. 

We also considered ‘reflection models’ of the RXTE data. 
These models showed evidence of a hardening of the RXTE 
spectra with decreasing X-ray flux. They further showed evi- 
dence of a positive correlation between the best-fit ionization 
parameter, f, and the observed flux. Especially the latter of 
these correlations, however, was dominated by the model fits of 
the faintest observation. The reflection models did not exhibit 
any evidence of a correlation between the photon index of the 
incident power law flux and the solid angle subtended by the 


reflector. 

Three of the RXTE observations were strictly simultaneous 
with 843 MHz and 8.3-9. 1 GHz radio observations. The most 
likely souice of the radio flux is synchrotron emission from 
an extended outflow with extent of 0( 10 7 GM/c 2 ). The cor- 
relation be:ween radio and X-ray emission on timescales of 7 
days or less (Hannikainnen et al. 1998) implies a strong cou- 
pling of tie inner disk accretion flow with this spatially ex- 
tended out low. Recent theoretical arguments have shown that 
one generally expects outflows in ADAF models (Blandford & 
Begelman 1998). Further simultaneous radio/X-ray observa- 
tions, preferrably with the addition of IR/optical monitoring to 
constrain the location of the synchrotron break and with the ad- 
dition of soft X-ray monitoring to constrain the accretion disk 
parameters , are required to test such models in detail. 
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We extracted data from the two solid state detectors (SISO, SIS 1 ) and the two gas detectors (GIS2, GIS3) onboard ASCA by using 
the standard ftools as described in the ASCA Data Reduction Guide (Day et al. 1998). The data extraction radius was limited by the 
fact that all the observations were in 1-CCD mode and that the source was typically placed close to the chip edge. We chose circular 
extraction regions with radii of as 4 arcmin for SISO, « 3 arcmin for SIS 1 (the maximum possible extraction radii for these detectors), 
and « 6 arcmin for GIS2 and GIS3. For observation 3 we excluded the central 40arcsec to avoid the possibility of photon pileup. 
We used the sisclean and gisclean tools (with default values) to remove hot and flickering pixels. We filtered the data with the strict 
cleaning criteria outlined in the ASCA Data Analysis Handbook; however, we took the larger value of 7 GeV/c for the rigidity. We 
rebinned the spectral files so that each energy bin contained a minimum of 20 photons. We retained SIS data in the 0.5 to 9 keV range 
and GIS data in the 1 .5 to 9 keV range. 

We accounted for the cross-calibration uncertainties of the three instruments by introducing a multiplicative constant for each 
detector in all of our fits. Relative to SISO, the SIS 1 detector normalization was always found to be within 2%, the GIS2 normalization 
was found to to be within 9%, and the GIS3 normalization was found to be within 15%. For any given observation, the normalization 
constants varied by ^ ±1% for different spectral fits. The background was measured from rectangular regions on the two edges of 
the chip farthest from the source (SIS data), or from annuli with inner radii > 6 arcmin (GIS data). These data were cleaned and 
filtered in the same manner as the source files. 

The resulting data files showed reasonable agreement between all four detectors. The most discrepant detector was SIS 1 , which 
also was the detector limited to the smallest extraction radius. This detector tended to show deviations from the other detectors for 
energies £ 9 keV, and from the SISO detector for energies « 0.5-1 keV. The detectors were mostly in mutual agreement for the lowest 
flux observations. It is likely that the agreement could be further improved for observations located closer to the center of the chips 
(thereby allowing larger extraction radii), and if low galactic latitude dark sky observations in 1-CCD mode were available to use as 
background. 


RXTE DATA EXTRACTION 

In §3 we present data from both pointed instruments on RXTE, the Proportional Counter Array, PCA, and the High Energy X-ray 
Timing Experiment, HEXTE. As we showed in §3, the large effective area of these instruments results in a data analysis approach 
that is dominated by the calibration uncertainty of these detectors (especially the PCA). In this appendix we summarize the major 
properties of both instruments and study their (inter-)calibration. All RXTE results obtained in this paper were obtained using 
the standard RXTE data analysis software, ftools version 4.1 (including the xte-patch 4.1 1 and the correct accounting of the time- 
dependence of the PCA response; Jahoda, 1998, priv. comm.), in our analysis. Spectral modeling was done using XSPEC, version 
10.00s (Amaud 1996). 

The PCA consists of five nearly identical co-aligned Xenon proportional counter units (PCUs) with a total effective area of about 
6500cm 2 . The instrument is sensitive in the energy range from 2keV to ~ 60keV (Jahoda et al. 1996). We only used data where 
all five PCUs were turned on. Background subtraction was done in the same manner as for our RXTE Cyg X-l observations (Dove 
et al. 1998). Specifically, a model using the rate of Very Large Events in the detector was used to estimate the background flux. 
The major uncertainty of this estimated background is due to activation of radioactivity in the detectors during SAA passages. Since 
this background component is present for about 30 minutes after the passage, we ignored data measured during these intervals. 



Fig. B8. — Ratio between the power-law fit of eq. (Bl) and the PCA Crab nebula data. The major deviations between the data and the model are 
due to the Xe edges at ~ 5.1 keV, and due to the Xe K edge at 34.6 keV. 
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Furthermore, data were not accumulated at times of high electron contamination. The electron contamination is measured by a 
certain ratio of veto rates in the detectors, the so-called “electron ratio” (Jahoda et al., 1998, in preparation). As recommended by the 
RXTE Guest Observer’s Facility (GOF), we excluded times during which the “electron ratio” was larger than 0.1 in at least one of 
the detectors. Note that the observed count rates from GX 339-4 are too high tc allow the use of the newer background model made 
available by the RXTE GOF in 1998 June. 

For spectral fitting, we limited the energy range of the PCA data ffom 3 to 30keV and used version 2.2.1 of the PCA response 
matrices. These matrices are newer than those used by us previously (Dove et al. 1998), and they are primarily characterized by 
assuming a higher instrumental resolution (Jahoda et al., 1998, in preparation). Due to the large PCA count rate of GX 339-4 
(~ 800 cps) our observations are dominated by the remaining uncertainties in the detector calibration, and not by Poisson errors. 
Therefore, a good understanding of these uncertainties is necessary. 

Since the Crab spectrum is commonly assumed to be a pure power-law (Toor & Seward 1974), at least over narrow energy ranges, 
the ratio between the fit to the Crab and the data can be used to deduce the systematic uncertainty associated with the detector 
calibration. We therefore extracted a public domain spectrum of the Crab nebula and pulsar measured with the PCA on 1997 April 1 . 
The Crab data were screened using the same criteria as those applied to our GX 339-4 data, except for that the “electron ratio” check 
was not applied since the background contributes only 0.6% to the total number of photons detected for the Crab. Modeling the 3 to 
30 keV Crab data with an absorbed power-law resulted in a best-fit photon spectrum of the form 

Aph = 13.3JT 2 187 exp [-2.54 x 10 21 cm* 2 <7bf(E)] ernes' 1 keV^ 1 (B1 ) 

where cr^E) is the energy dependent bound free absorption cross section for material of cosmic abundances as given by Morrison 
& McCammon (1983), and where E is the photon-energy measured in keV. For this fit to the Crab data, x 2 /dof = 168^56. In Fig. B8 
we display the ratio between the best-fit to the Crab and the data. Using this ratio plot, we deduced the systematic uncertainties of 
the detector (Table B9). Adding them in quadrature to the Poisson errors of our data, the x 2 of the above Crab fit was reduced to 
Xred = 0*56. Note that the internal consistency of the detector calibration of the PCA appears to be on the 1% level, i.e., comparable 
to that obtained for previously flown instruments like the Ginga LAC (Turner et al. 1989) or the ASCA GIS (Makishima et al. 1996), 
even though the much larger effective area of the PCA makes it necessary to include many “dirt effects” into the detector model. 

We caution that the statistical uncertainties of the fit parameters presented in this work were derived using the above systematic un- 
certainties. It is questionable, therefore, whether the approach of Lampton, Margon & Bowyer (1976) to determine the uncertainties 
from the x 2 -contours can really be used since this approach makes use of the assumption of Poisson-type errors. The uncertainties 
given in this work should be taken with these caveats in mind, especially for hose fits where the x^a values are very small (i.e., 
£ 0.5), and thus they should not be construed literally as ‘90% uncertainties’. 

The HEXTE consists of two clusters of four Nal/CsI-phoswich scintiliatior counters, sensitive from 15 to 250 keV. A full de- 
scription of the instrument is given by Rothschild et al. (1998). Source-background rocking of the two clusters provides a direct 
measurement of the HEXTE background with measured long term systematic uncertainties of < 1% (Rothschild et al. 1998). Al- 
though no other strong sources in the field around GX 339-^1 are known (Covault, Grindlay & Manandhar 1992; Bouchet et al. 1993; 
Trudolyubov et al. 1998), we extracted individual background spectra for both HEXTE cluster background positions to check for 
contamination of the spectrum ffom weak background sources. In all cases the background spectra differed by £ 1 cps. Thus, we 
used the added background spectra ffom both cluster positions in our data analysis. We used the standard response matrices dated 
1997 March 20 and used data measured between 17 and 1 lOkeV. An analysis of the detector calibration similar to that performed for 
the PCA reveals that the HEXTE calibration is good on a level comparable to th ; PCA. Due to the much smaller effective area of the 
detector and due to the smaller flux from the source at higher energies, howeve*, the HEXTE spectrum is completely dominated by 
the Poisson error of the data. Therefore we did not attempt to take the systematic calibration uncertainty into account. To improve 
the statistics of individual energy bins, we rebinned the raw (as 1 keV wide) er ergy channels by a factor of 2 for the energy range 
from 30 to 5 1 keV, and by a factor of 3 above that. 

When modeling the spectrum of GX 339-4 from 3 to 1 lOkeV as measure! by the PCA and the HEXTE, the intercalibration 
between the instruments is of some concern. Our experience with previous daa and the response matrix described by Dove et al. 
(1998) indicated that the flux calibration of the HEXTE with respect to the PCA was off by about 25%, i.e,, the derived HEXTE 
fluxes were ~ 75% of those found with the PCA. This deviation is mainly due to a slight misalignment of the HEXTE collimators 
which has not yet been taken into account in the HEXTE response matrix (Heintll 1998, priv. comm.). Using the new PCA response, 
the flux ratio now appears to be larger and was found to be 62 — 69% in our data. Extracting spectra with internal software used by 
the HEXTE hardware team produced spectra identical to those found using ou: extraction procedure. Therefore we do not believe 
this change in the flux calibration to be due to errors in the deadtime correction. To take this offset in the effective areas into account 
we modeled the spectra using a multiplicative constant which was set to unity ) or the PCA, and which was a fit parameter for both 
HEXTE clusters. Thus, all fluxes given below were measured with respect tc the PCA. The maximum deviation of the HEXTE 
clusters relative to each other was found to be less than 8%. 


Table B9: Systematic uncertainties of version 2.2.1 of the PCA response matrix as obta ned from a power-law fit to the Crab spectrum (Fig. B8). 

Channel 1 0-15 16-39 41^-57 58-128 

Channel-Energy [keV] 0-7 7-16 l<-25 25- 

Uncertainty [%] 1 0.5 2 5 


PCA channels assuming the standard2 channel binning. 
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Apart from the flux calibration, however, even more crucial for our analysis is the question of how well the inferred spectral shapes 
agree for the two instruments. Our Crab fits show that the PCA results in a photon index of T = 2. 1 87 (Eq. Bl), while our HEXTE 
fits gave T = 2.053. The generally accepted value for the Crab photon index in the 1-100 keV range is T = 2. 10 ± 0.03. This value 
was adopted by Toor & Seward (1974) in their analysis of 28 different rocket flight measurements. There are indications that the 
spectrum softens to T ~ 2.5 above 150keV (Jung 1989). Although the absolute uncertainty of the Crab flux in the 2-100keV range 
has been estimated as large as 2%, and even larger below 2keV (Nprgaard-Nielsen et al. 1994), the AT = 0.134 deviation between 
the PCA and the HEXTE best-fit Crab photon index is still very worrisome, and it is currently being studied by both the PCA and 
HEXTE instrument teams (Jahoda 1998, priv. comm.). 
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